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It is known from the LighthilLs exact solution of the incompressible inverse problem 
that in the inverse design problem the surface pressure distribution and the free stream 
speed can not both be prescribed independently. This implies the existence of a constraint 
(regularity condition) on the prescribed pressure distribution. The same constraint exists 
at compressible speeds. In this paper, a well-posed inverse design method for transonic 
airfoil is presented. In the method, the target pressure distribution contains a free 
parameter that is adjusted during the computation to satisfy the regularity condiuon de- 
rived in this paper. A few design results are presented here in order to demonstrate the 

capability of the method. 


L 


INTRODUCTION 

Recently, a number of design methods have been developed and used for the design 
of transonic airfoils and wings. Slooff 1 reviewed these methods and divided them into 
three major categories: indirect, inverse, and aerodynamic optimization. Indirect meth- 
ods are characterized by the fact that the designer has no control over either the 
aerodynamic quantities or the geometry. The hodograph and fictitious gas methods are 
in this category. In inverse methods, the classical inverse problem of aerodynamics is 
solved. The designer specifies an arbitrary pressure distribution on an airfoil or wing, 
while the geometry of the airfoil or the wing that realizes the given pressure distribution 
is determined as the result of the solution. Aerodynamic optimization methods are those 
in which a nonlinear optimization algorithm is linked with a flow analysis code to mini- 
mize or maximize some aerodynamic object functions such as the lift-to-drag ratio. 

The conventional inverse design methods are the most used in the industry app ica 
tion. The currently existing inverse methods for transonic airfoil design can be su i- 
vided into two categories: (a) methods utilizing Dirichlet-type boundary conditions de- 
rived from the target pressure distribution; (b) methods utilizing Neumann type oun 
ary conditions in combination with some geometry correction procedure (residual- cor- 
rection method). . „„ 

In fact in inverse problems, both a Dirichlet- and a Neumann-type boundary con- 

dition must be satisfied on the airfoil contour to be determined. This gives rise to a 
nonlinear problem with unknown boundary to be solved iteratively. n 1 e irs 
approach, the required target pressure distribution is imposed on an initial airfoil as a 
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Dirichlet boundary condition, while airfoil geometry corrections are derived by inte- 
grating (either explicitly, or in some implicit manner) the transpiration mass flow over 
the initial airfoil. The Neumann boundary condition is satisfied at the end of the 
iterations. 

In the second approach, the pressure distribution on an initial airfoil is determined 
by the use of an analysis code (Neumann boundary condition), and the residuals (i.e. dif- 
ferences between target pressure distribution and pressure distribution on the airfoil of 
the current iteration) are transformed to airfoil geometry corrections using some relative- 
ly simple approximate inverse methods. Here, the Dirichlet boundary condition is satis- 
fied at the end of the iterations. 

The main advantage of the iterative Dirichlet-type method is rapid convergence, 
provided closure and regularity conditions are satisfied. This kind of methods has been 
developed by Volpe and Melnik, 2 Carlson 3 and Tranen 4 for airfoils, Henne 5 for wing, 
and Shankar 6 for wing / body designs. But in most of these methods, the regularity con- 
dition is obviously not taken into account. 

The main advantage of the residual-correction method is its simplicity. Only a small 
investment is required, because any existing transonic analysis code can be used without 
modification. Efforts can be concentrated on coding a simple approximate inverse rou- 
tine base on a suitable approach, and coupling this to the analysis code. Another advan- 
tage of this kind of methods is that the analysis code can be easily replaced with more 
advanced codes when they become available. This approach was used by 
Davis, 7 McFadden, 8 Fray el al 9 Greff and Mantel. 10 Takanashi" has developed a 
three-dimensional transonic wing design method based on this approach. The geometry 
correction problem is formulated in a three-dimensional transonic small— disturbance in- 
tegral equation form and is numerically solved. 

It was demonstrated by Lighthill 12 that a unique and correct solution to the inverse 
problem of 2— D, incompressible flow generally does not exist unless the prescribed speed 
distribution satisfies a certain integral constraint (regularity condition) arising from the 
requirement that the speed in the free stream be equal to one (or any other specified val- 
ue). Thus, in order to assure that a solution to the inverse problem exists, some freedom 
must be permitted in the prescribed pressure distribution to allow the regularity condi- 
tion to be satisfied. This can be accomplished by introducing a parameter into prescribed 
pressure distribution. In transonic flow the similar constraints have never been properly 
formulated. Volpe and Melnik 2 proposed an inverse design method for transonic flow 
which was aimed at providing a treatment of such kind of constraints on the target speed 
distribution numerically. 

In this paper, we formulate a well-posed inverse design method for transonic airfoil. 
This method is a residual-correction type approach, in which the procedure to determine 
the geometry correction is similar to the procedure used by Takanashi for his 
three-dimensional wing design method. An integral constraint (regularity condition) is 
derived and used in the method. Thus, a parameter is introduced into the target pressure 
distribution, and this parameter is determined as part of the solution according to the 
regularity condition. 
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FORMULATION OF THE INVERSE PROBLEM 
The full potential equation can be written in terms of the perturbation velocity po- 
tential as 

' 2 4 - m — k 

dx 

Where M^ is the free stream Mach number, K^M,*, ) is a transonic similarity 
parameter, y is the ratio of specific heats, (x,z) is a Cartesian coordinate system, and H 
represents all the higher order terms. 

the tangency condition on the airfoil may be written as 

V; (x,±0) = f' ± (x) + Q (2) 

where f (3c) and f _ (3c) are the equations of the upper and lower airfoil surfaces 

respectively, and Q represents all the higher order terms. 

the pressure coefficients on the airfoil surface are expressed as 

C (x)= -2 qAx,±Q)+S ( 3 ) 

where S represents all the higher order terms. 

In a residual-correction method, the solution <p(x,z) of Eq. (1) for an initial 
airfoil f (x) has been obtained by means of an existing analysis code, the objective here 

is to determine the amount of the geometry correction A? ± (x) corresponding to the 
pressure difference A C p± (x) between the specified and calculated pressure. If a small 
perturbation A (p(x,z) is further introduced, we can obtain the potential equations for 
A <p(x,z) according to the transonic small-disturbance theory 


A© + A© = 


dx 


+ A< iP x ) 2 -^x] 


A <p t (x, + 0) = Af ' ± (x) 

n 2 

AC p± (x)= — 2-g A(p x (x, ± 0) 


( 4 ) 

( 5 ) 

( 6 ) 


( 7 ) 


here /? = — M 2 , and the new variables have been introduced as 

x = x, z = Pz, <p(x,z)=(K/ P 2 )v(xj), f ± (x)=(K/ P 2 )l ± (F) 

By applying Green's theory to Eq.(4) and introducing a decay function similar to 
that used by N<t»rstrud 13 , we can get the integral equations 

A u (x) = -f ¥ (x,0;£,0)Aw (flrff + (?,(*) 
s nJ 0 * 

- - f [I (x£, + 0)G(f, + 0) + I (x,^ - 0)(7(x, - 0m 
71 J o ' 

r j Am ({) 


( 8 ) 


1 
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-- ['[/ (x^, + 0)G(^ + 0)-/ a (x^,-0)G(x,-0)]^ (9) 

7T J o 

here 

Au s (x) = A<p x {x, + 0) + A<p x (x, — 0) (10) 

Au o (x) = Acp x {x, + 0) — A(p x (x, — 0) (11) 

Ah’^jc) = A(p z (x, + 0) — Acp^x, — 0) (12) 

Aw a (x) = A(p z (x, + 0) + A(p z {x, — 0) (13) 

¥(*,*;{, 0 = ln[(x - O' +(z~ Q 2 ] 2 04) 

G(x,z) = ^ (<p x + A(p x f - <p\ J (15) 

I (x,£,±) = AxM,0exp[-2R ± «m 06) 

/.(*,{, + ) = | X V {z (x,0;£,O«/>[ - 2i? ± 07) 

-R ± (*) = |f'i (x)/<P x (x, + 0)| (18) 


For convenience, the correction function Af ± (x) is split into symmetric Af s (x) and 
antisymmetric Af a (x) parts 

Af^(x) = Af + (x) — Af _ (x) (19) 

Af a (x) = Af + (x) + Af _ (x) 

Since Af a (x) = Aw a (x), the antisymmetric part can be determined by direct evaluation 
of the right-hand side of Eq. (9). On the other hand, since Af s '(x) = Aw s (x). the symmetric 
part must be solved implicitly. Consequently, the correction Af ± (x) is obtained by inte- 
grating Aw s (x) and Aw a (x) with respect to x. 

CONSTRAINTS FOR INVERSE PROBLEM 

In the incompressible flow, it was demonstrated by Lighthill, 12 using conformal 
mapping method, that the geometry of an airfoil for a given speed distribution can be de- 
termined only if the prescribed speed distribution satisfied the following three integral 
constraints 

log \— — | do = 0 (20) 

J o Q & 

f 2 * d n 

log | 1 COSO) do) = 0 (21) 

o d a, 

log\— | sirno do = 0 (22) 

J 0 d oo 

here q 0 is the prescribed speed distribution on the airfoil surface, q* is the speed at infin- 
ity, o is the polar angle in the transformed plane. The first constraint known as regulari- 
ty condition is a consequence of the fact that the speed at infinity is q^. Eqs. (21) and 
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^ (22) together express that the airfoil is a closed contour and the angle of incidence is zero 
(or any other specified value). 

In order to formulate the regularity condition for transonic flow, differentiating 
both sides of Eq(4) with respect to x, we have 


V V(x,z) = 


dx 


here 


V(x,z) = Acp x (x,z) 
with a Dirichlet— type boundary condition 

VJx,±0)=V o (x)= 

and a constraint at infinity 

A(p x (x,z)^0 ; (x 3 

Now we use the transformation 

1 


T(x) =x + iz = -[ x + 


9*H( 


> 

(23) 


(24) 

AC, t (x) 

(25) 

2 


- Z —-+ oc) 

P 

(26) 

iio e' to \, 1 

" + r ) + 2 

(27) 


where x = re‘ w is the complex variable in the transformed plane. The entire plane in the 
physical plane is mapped onto the outside of the unit circle, the chord of the airfoil (z = 0, 
O^x^l) is corresponds to the unit circle (r=l, 0^<u^7i) on the transformed plane. 
Substituting Eq . (27) into Eq . (23) yields 

V 2 V(r,oj) = L / \ T'(x)\ 


(28) 


where 


\ d 2 G i 

f dr' 

\8 G 

dr 

dG 

d co 

dG 

)drdco + 1 

v dx - 


+ dx 2 

dr 

+ Sx 1 

dco 


dx J do 2 Vdxi 

The boundary conditions in the transformed plane are 

V(r,co) | r _, = V 0 (co) 

V ->0 ; (r-*oo) 

By applying Green's theory to Eq.(28), an integral expression can be obtained as 

rj-lr 2 * V 0 (a>) 

V{r , ,<u j ) — J 0 r t — 2r ^cosico - co^ + 1 


(29) 

(30) 


+ 2n 


In 


log 


r j — 2 rr ( cos(co — co ] ) + r 
Lr 2 r — 2 rr cos( co — co^ + 1 J I T\x)\ 


rdrda 


If r — » oo, the above equation becomes 

In physical plane, Eq. (32) can be written as 

f r a m (y -i- O'! _ Am ( r — 0)lru dx = J J logr 


g 2 G(x,z) 
a * 


dxdz 


(31) 


(32) 


(33) 
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Using an integration by parts and introducing the decay function, we can obtain the 
final form of the regularity condition 

1 r 1 

A u a (x)co x dx= [I'(x, + Q)G(x, + Q) + I y (x,-0)G(x,-0)]dx (34) 


here 


/ Ax, ± 0) = 


3 rr 


-exp[ — 2R ± (x)z]dz 


(35) 


o r 

The above discussion indicates that the prescribed pressure distribution should con- 
tain an adjustable parameter to guarantee that the regularity condition is satisfied. Thus 
the surface pressure distribution is to be prescribed in the form 

C ps± (x) = F ± (<r,x) (36) 

where cr is a parameter that is found as part of the solution. For convenience, the follow- 
ing form of C DS is used in this paper: 

c ps± (x) = cF ± (x) (37) 

In order to assure that the resulting airfoil has a specified trailing edge gap, the fol- 
lowing closure condition must be satisfied: 

• i 

Avv j (x)i/x = 0 (38) 

o 1 

This closure condition can assure that the trailing edge thickness of the current 
airfoil is always kept equal to that of the initial airfoil. 

Although the regularity condition in closed form (34) is obtained from the simple 
transonic small disturbance theory with N<)>rstrud assumption, its practical utility will be 
shown by numerical examples in the following section. 

NUMERICAL AND COMPUTATIONAL ASPECTS 

In order to discrete the integrals appearing in Eqs. (16) and (17), the range of 
integration with respect to ( is divided into subintervals. Assuming that R ± (£)( ‘ s con- 
stant on each subinterval, Eqs. (16) and (17) can be expressed as 


l 




(«.]J 

c4 AC - 

f C. 

CM 

(39) 

KXJ_ 

c„-X 


(40) 


The range of integration with respect to x is also divided into subintervals, and on 
each of subinterval A u } (x), A u a (x), G(x,± 0), and Aw a (x) are assumed to be con- 
stants, while A w s (x) is assumed to vary linearly, the final expressions of Eqs. (8) and (9) 
in discretized form are as follows: ^ 

Au a (x.)= Y. d ik Aw s (x k _\ ) + G 3 (x.)- Z [v’ t G(x t , + 0) + v rt G(x t ,-0) 

4-1 2 1 


(41) 


L 
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A w (x )= £/i‘ t Au a (x t )- X v ° t G(x t , + 0) - v. k G(x k , 0)j 

* 4 4-1 

The discretized form of Eq. (34) and (38) can be expressed as 
XAu a (x. )(<u. + i ) = Z v'G(x., + 0) + v'G(x.,-0)J 


(42) 

(43) 

(44) 


I [Aw,(TJ ) + Aw j (x i t )](* i+ l ) =0 

I- 1 2 

where the coefficients 

s a s *3 a r 

^ i* ’ ^ ik * V i* * V ik 9 V ik’ V ik 9 i* i 

are the integral expressions on each of the subinterval. The full expressions are omitted 

here. , 

In order to improve the convergence, the following modifications have been taken m 

this paper: 

1) . A Riegels type of leading edge correction is taken in the method. The purpose ot 
such correction is to remove the singularity at the leading edge of the round-nosed 

airfoil. . . 

2) . In order to increase the ability to deal with the shock, an artificial viscosity term 

is added to the integral equation method. 

3) . A Smoothing-relaxation procedure is proposed and used in this paper. 

The inverse problem can be solved by the iteration process as follows: 

1) The flowfield is solved for an initial airfoil f ± (x) by a direct analysis code. From 

the calculated pressure distribution C p± (x), and the target pressure distribution 
C DS ±(x), the residual AC p± (x) = ffC ps± (x)-C p± (x) can be obtained. 

2) . The adjustable parameter a is determined from the regularity condition (43). 

3) . The geometric correction Af ± (x) is determined by solving the equations 
(41)-(42). Thus a new geometry is obtained from the following Smoothing-relaxation 

procedure: . \ q 

f*/ 1 !*,)- r* ± (*,) + <[ AC l (*.) + 0 - 5 ( Af *f ,(x '-' ) + AI * 1 (x <*> ) )J (45) 

where <5 is a relaxation factor. 

The same process is repeated until the calculated pressure distribution agrees wit 
the prescribed one. 

Several test cases are presented to show the validity and applicability of the proce- 
dure. A nonisentropic potential solver for 2-D transonic flow' 4 is used as the analysis 


In the first example, the target pressure distribution is taken from the result of 
RAE2822 airfoil, the free stream Mach number is 0.73, and the angle of attack is 
ot = 2 05 ° . The initial airfoil is NACA0012 airfoil, the initial angle of attack is a -0.0. 
After 12 design cycles, the RAE2822 airfoil is already recovered. Indeed the designed 
airfoil is rotated in a clockwise direction by 2.05 ° with respect to the original airfoil. 
This is because the initial angle of attack is chosen to be 0 ° , the angle between the free 
stream and the X-Axi is always 0 ° during iterations. Fig.l shows the target and initial j 
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pressure distribution. Fig. 2 shows the convergence procedure of the pressure 
distribution. 

The second example is chosen to show the function of the regularity condition. In 
Ref. [15], Strand gave a velocity distribution which did not satisfy the LighthilTs three in- 
tegral constraints in incompressible flow, hence there is no airfoil corresponding to this 
velocity distribution. Strand modified this velocity distribution according to the 
Lighthill's constraints, and designed the airfoil corresponding to the modified velocity 
distribution, the results is showed in Fig. 3. Now, we calculate this example using the 
present method. For convenience, only the velocity distribution on the lower surface is 
multiplied by the free parameter cr. After 7 design cycles, the converged solution is ob- 
tained, the velocity distributions are given in Fig. 4. In this case, the output velocity dis- 
tribution on the lower surface is not consistent with input velocity distribution, this is be- 
cause the free parameter a is not equal to 1, but equal to 1.07. The results show that the 
method can adjust the improperly input velocity distribution to the acceptable velocity 
distribution automatically, and design the airfoil corresponding to the modified velocity 
distribution. Comparing the figures the present result is closer with input data than 
Strand's one. On the other hand, if the regularity condition is not included in the 
method, the design procedure will not converge for this example. 

CONCLUSIONS 

A regularity condition in closed form for transonic flow is presented in this paper, 
and a well posed inverse design method for transonic airfoil is formulated. The results 
show that the method is a reliable and efficient method for the design of airfoil at 
transonic speeds. When the target pressure distribution is not properly given, the code 
can adjust the target pressure distribution automatically and design the airfoil corre- 
sponding to the modified pressure distribution. 
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Fig. 2 Convergence history of the pressure distribution (exatple 1) 


Cp A -A Initial 




551 

Third International Conference on Inverse Design Concepts and Optimization in Engineering Sciences 
nCTDES-IIIi. Editor: G.S. Dulikravich. Washington P.C.. October 2 Ii2L 122L 


Turbulent-flow 



Distance along surface, leasured fron 
leading-edge stagnation point 
Fig. 3 The input and acceptable velocity 
distribution given by Strand 
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Fig. 4 The input and designed velocity 
distribution (example 2) 
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